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Abstract 



^\j , Recently, Belle and BaBar Collaborations observed surprising suppression in the 

endpoint J/ip spectrum, which stimulates us to examine the endpoint behaviors of the 

e~^e~ — > J/ipgg production. We calculate the J/ip momentum and angular distributions 

■r^ij- . for this process within the framework of the soft-collinear effective theory (SCET). 

^^ I The decreasing spectrum in the endpoint region is obtained by summing the Sudakov 

,^ ' logarithms. We also find a large discrepancy between the NRQCD and SCET spectrum 

'j^. in the endpoint region even before the large logarithms are summed, which is probably 

Qh! due to the fact that only the scalar structure of the two-gluon system is picked out in 

^ ' the leading power expansion. A comparison with the process T — > ^gg is made. 

> 
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1 Introduction 

Heavy quarkonium system plays an important role in the development of quantum cliromo- 
dynamics (QCD). The scale of the heavy quark mass guarantees the applicability of pertur- 
bative QCD, meanwhile the nonperturbative physics presents itself through hadronization 
effects. In the past few years, one of developments in heavy quarkonium physics, called 
nonrelativistic quantum chromodynamics (NRQCD) PP which generalizes and improves the 
conventional color-singlet model (CSM), has provided a successful explanation of the sur- 
prising excesses of J/V' and ip' productions at the Tevatron |2] by introducing color octet 
contributions. 

NRQCD factorization should be further examined in other collider facilities, in particular, 
e^e~ colliders which provide a clean testing ground. SLAC and KEK e"^e~ B factories are 
now running at or below the T(4s) resonance. At this energy, it was expected in NRQCD 
that the inclusive J /ipgg process should be dominant |S1I1] and in the upper endpoint region 
of the J /ip momentum spectrum, there may exist a sharp peak as a clean signal of the 
color-octet ccg contribution j3]. 

Recently, BaBar [H] and Belle [3 |H] Collaborations published their measurements on 
prompt J /ip productions in e^e~ collision at center-of-mass (cm.) energy -y/i = 10.58 GeV. 
It is really surprising to observe that, according to Belle's data [7||H], it is the J /ipcc process 
that dominates the inclusive J/ip production at B factories 

a(e+e' -^ J li)cc)lo{e^e- -^ J/ipX) = 0.67 ± 0.12, (1) 

while the momentum distribution of the inclusive J/ip production shows a suppression, 
instead of an (expected) enhancement, in the upper endpoint region. 

For the unexpected J/ipcc dominance, it is argued in Ref. j3] that a large renormalization 
K factor might be the answer. Recent investigation [9^ also reveals that the color-octet contri- 
bution to J/ip spectrum can be broadened significantly by the large perturbative corrections 
and enhanced nonperturbative effects so as not to conflict with the surprising suppression in 
the endpoint region observed by BaBar and Belle. However a leading-order NRQCD calcu- 
lation shows that, in the endpoint region, the color-singlet J/ipgg contribution is not small 
at all, which seems to be still in contradiction with the experimental observations. In this 
work, we are stimulated to investigate the endpoint behaviors of the e+e~ — >• 7* -^ J/'Pgg 
production. 

We note that, at the amplitude level, 7* —>■ J/ipgg is very similar to the decay T —>■ '~fgg. It 
has been known several years ago that, at the endpoint of the photon spectrum in radiative T 
decays, NRQCD is not applicable due to the breakdown of both the perturbative expansion 
and the operator product expansion (OPE) JUl- The same arguments should also apply 
for the case of the J/tp production. This is because NRQCD only contains soft degrees of 
freedom at low energy, but at the endpoint of the photon and/or J/ip spectrum, the gluon 
jet should be almost collinear. To fix this problem, Fleming et al. proposed a combination of 
NRQCD for the heavy degrees of freedom and the soft-collinear effective theory (SCET) JT] 
for the light degrees of freedom. With this method, the radiative T decays were investigated 



in a series of papers |T21 El El which show an improved agreement with the CLEO data 
[T3] . Lately the same method was apphed to the color-octet contribution to the inclusive J/ifj 
production e^e~ -^ J /tp + X jH]. By the use of the resummation of Sudakov logarithms and 
the nonperturbative shape functions, the color-octet J /ip spectrum, which is a sharp peak at 
maximal energy in leading order calculations, could be significantly broadened and shifted 
to lower energies. This therefore would resolve the discrepancy between the color-octet 
J /ip production and the experimental observations. According to the spirit of the Sudakov 
suppression in the endpoint region, the authors in Ref. |3] adopted a phenomenological 
approach to obtain an appropriate endpoint spectrum for the J/ipgg process instead of 
performing a complete calculation in SCET. 

In this paper, we shall follow the same way of Ref. ^31^]; namely SCET combined with 
NRQCD, to examine the endpoint behavior of the color-singlet J/ipgg mechanism. 

2 Leading order SCET calculation 

Several scales are involved in this process: the center-of-mass energy ^/s, the J/tp mass 
M^, and the non-perturbative QCD scale Aqcd- In this paper we will only consider the case 
where the ratio M^j ^ is kept finite in the limit of infinite y^. In this point of view, J jip can 
be taken as a heavy particle. In the kinematic endpoint region of J jip spectra, the failure of 
NRQCD factorization and the relevance of SCET has been explained clearly in Refs. Plll3j. 



In brief, the hadronic jet recoiling against J j^tp is not highly virtual, rrtx ~ \j ^fs^QCD-, 
compared with its large momentum of order -y/s. This results in the OPE breaking down, 
and therefore a new effective theory, the so-called SCET, is developed by including collinear 
degrees of freedom. 

In SCET, it is convenient to write a momentum in light-cone coordinates. Working in 
the e^e~ cm. frame, we define the incoming electron and positron moving along light-cone 
directions n^ = (1, 0, 0, —1) and n^ = (1, 0, 0, 1). The produced J/ip meson is chosen to move 
in the x — z plane with momentum Pjt = Mv^ = {E, \P\ sin6, 0, \P\ cos6) (M is J/ip mass), 
and hence the light-cone vectors for two gluons can be defined as n'^ = (1, — sin^, 0, — cos6') 
and n'^ = (1, sin6', 0, cos^). Throughout this paper, we adopt a dimensionless variable z = 
\P^\/P^°'^, where P^""^ denotes the maximum value of the J/ip momentum, namely PJ^"^ = 
■\/s{l — r)/2 ~ 4.9 GeV. Here r = M'^/s ~ 0.08. The J/ip velocity v can be expressed as 

v'' = (vo, ItTlsin^, 0, Iwlcos^) = (\ r ~ ^^ z^ + 1 , —^zsmO, 0, —^zcosO). (2) 

V 4r 2y/r 2y/r 

For the process e~^e~ — ;► 7* — * J/ipX, the hadronic jet has the momentum px = l^ — Mv^ — 
k^^ where l^ = {^/s, 0, 0, 0) is the momentum of the virtual photon and k^^ is the residual 
momentum of the cc pair within J/ip. In the endpoint region, since the hadronic jet is 
collinear along the light-cone direction n'^, we can write px ~ V^(l; A^, A) in the n — n light- 
cone coordinate. When E™"^ — E^ ~ Aqcd, Px is of order 2A/i(-E'™"^ — E^) ~ 2^/sAQCD 
which implies NRQCD factorization breaks down in this kinematic region. Therefore SCET 



becomes relevant in the endpoint region 1 — 2; ~ Kqcd/M ~ f^, and correspondingly the 
expansion parameter A is of order a/1 — z in this process. 

Before going into details, it is helpful to notice the similarity between e+e^ — * 7* — >^ 
J/ipgg and T — >• ^gg. In fact the cross section of J/ipgg production can be related to the 
"decay width" of the transversely polarized virtual photon 7* as follows [12] 

dcr(e+e- -^ J/i^gg) = A7ias-^/^dT{-f* -^ J/i^gg) , (3) 

where the polarization vector of the virtual photon satisfies the following equation 

e^e*'^ = ~g^^ + ^-^-^^-^^ ^ -^ . (4) 

It is then clear that, at the amplitude level, the effective operator of 7* -^ J/'4'gg should be 
formally the same as that of T — i> 'jgg. Therefore the proof of SCET factorization for the 
former process is almost the same as that of the latter one, which has been elaborated in 
Ref . ^3] ■ All of our following calculations will be in parallel with those for T — > 'jgg in Ref . 

m 

To proceed, we shall first match from NRQCD onto SCET. Considering the gauge and 
reparametrization invariance, the leading SCET color-singlet ^Si operator is given by |13| ^ 

= ^lA . ^'x^,Tr{BlT%';\^^^, -n ■ vV^)B^,} , (5) 

where ipp and X-p ^-i'^ the heavy quark and antiquark fields from NRQCD, and B_j_ is 
the leading piece of the coUinear-gauge invariant gluon field strength jTHj. The operator 
V{V^) projects out the large light-cone momentum components of the collinear fields to 
the right (left). The second line of the above equation is obtained by using the identity 
Bfn ■ v{V + V^)B^ = -M(l - r)/rBlB^ and the definition P_ = P - pt. prom the 
matching shown in Fig. 1, we obtain the coefficient 

j.ii-'Si)f M{l -r) _ Ag^^eec r j_ 1-r _ 



r 



v/6 M(l-r)^"'''^'"' 2r 



where (7^^ = gap — ijiafip + npna)/^, n ■ q^ = n ■ q — n ■ q' . q and q' are the momenta of two 
gluons. 

According to the optical theorem, the J /ip momentum spectrum and angular distribution 
can be expressed as 

^ikL^i^ . (^r?^ i„r(.. 0) (7) 

dzdcO^e StT^^JM^ + (P"'^^)222 



^Here 'leading' means that the leading-order power expansion in terms of the small parameter A in SCET. 




+ crossed diagrams 

Figure 1: Matching the amphtude for 7* -^ J/'ipgg process in QCD and SCET. 

where the forward scattering amphtude is 

T{z,e) = -z fd'xe-^'^Y.{0\Jlix)\J/^ + X){J/^ + X\J,{0)\0)gZ. (8) 

In SCET, the following factorization formula can be proved in the endpoint region 

X J dk+Sik+, /i)Im^(A;+ + v^ - P;^'^";^ - ^ M^ + (P^»-)22^ /i) , (9) 

where H, S and J^ are the hard function, ultrasoft function and jet function, respectively. 
In order to obtain the above formula, we match the QCD current J^ in Eq. (jHI) to the leading 
SCET color-singlet operator Eq. © 

J,(x) = 5:e-(^^^-^W2))-.,r(,i^;ji)(^)^t^A ■ ^^X-pTr{P^5^,^_P^} , (10) 

where the operator V in the phase factor will sum the label momentum of the two collinear 
fields B± and thus can be replaced by —^5(1 — r). The matching coefficient '^a^St^^ i^) i^ 
given in Eq. (jH)). Since collinear fields in SCET are decoupled from ultrasoft gluons by field 
redefinition [llj and J ftp meson has no collinear freedom, the forward scattering amplitude 
in Eq. (jHl) can then be factorized by separating the heavy quark fields into ultrasoft functions 
and the collinear gluon fields into jet functions. Specifically, the jet function is defined from 
the vacuum matrix element of the collinear fields, which is exactly the same as that of the 
color-singlet radiative T decay [12] 

(0|TTr[i?f5^,P_Pf](x)Tr[Pf'5^,,P_i?f](0)|0) 

- liaT'gf + gfgi')6.,.' J ^.^-''-^Uk\ /i) , (n) 



where Bj_ is the redefinition of the colhnear field to decouple from the ultrasoft gluons. 
To calculate the jet function, one may directly evaluate the vacuum matrix element of the 
collinear fields, which is the left-hand side of Eq. (|TT| . Actually the jet function which is 
independent of the heavy quark fields, should be the same for both T — >■ 'jgg and 7* —>■ J/ipgg 
processes, so it can be obtained directly from Ref . [13] . For our purpose, only the imaginary 
part of the jet function is relevant, and at the lowest order in a^, it is 

Im^(fc+, /i) = ^0(^+) l\ d^S^^^s(i-r)i ■ (12) 

Following Ref. ^3, the ultrasoft function for this process can be written as 

= (0|xlpO-i^pa;);a^(5(m-(9-A;+)V'^o-iX-p|0) , (13) 

while the leading order hard function is computed as 



3 \^/6M{l 
327?^ / AasCecT 



2 



Fiz,e), (14) 



3 \^/6M{l-r], 
where the explicit expression for F{z, 6) is 

F(z, 6) =2- sin^ e + ^[(1 + r)v^ - (1 - r)\v\f . (15) 

Here v is the J /ip velocity given in Eq. (j2)). 

With these functions in hand, we obtain the explicit form for the imaginary part of the 
forward scattering amplitude in Eq. as follows 

1 (Stt / A(y pp T \ ± 



3M \^/6M{l- 



r\ 



W,|/?(0)p / 4a,ee,r 



2 



where -R(O) denotes the radial wave function of J/ip at the origin. Using the above equation, 
we arrive at the differential cross section in the tree-level SCET calculation ^ 

dzdcose ^^ ^ ^' 9s M3 /4^ + n - r)2^2 ^ ' ^ ^ ^ 



^'The tree- level SCET calculation' here is referred to the leading power calculation in SCET before the 
resumniation over large logarithms. 



Since OPE breaks down and large logarithms arise as z approaches to 1, resummation 
over large logarithms is indispensable before comparing to the experimental observations. 
In SCET, these logarithms can be resummed using renormalization group equations (RGE). 
To do this, one has to first calculate the anomalous dimension of the effective operator (Eq. 
(0)). Fortunately, this effective operator is formally the same as that which appears in the 
color-singlet radiative T decays, therefore we can directly read the anomalous dimension 
from Ref. JH] 



7W4{C. 



ii + (.= + (1 - .n f ^i^ + i^^^i^' 



With this anomalous dimension, one can then resum the large logarithms using RGE from 
the matching (hard) scale to the coUinear scale. The coUinear scale should roughly be the 



invariant mass of the jet, namely \ic{z) = J2^/s{E'!^°■^ — E^{z)). However there is no obvious 
clue what the matching scale should be. In NRQCD calculations, this scale is often chosen 
at quarkonium mass M, but in SCET it is found that, at least the hard scale for color-octet 
J/ip production should be about —n ■ vP = M(l — r)/r [0], according to the logarithm 
that appears in the anomalous dimension calculations. As we know that, there is no large 
logarithm far from the endpoint region, which means that the coUinear scale, which is of the 
order of ^/s for small z, should be comparable to the hard scale. Therefore in the following, 
we will naively choose the hard scale as fih = a/s(1 — i"), which is simply the large light-cone 
component of the gluon jet momentum. Finally, the result for the resummed differential 
cross section is 

daresnm r^r u ^ , 32(a,(/i;,)«e,)2Ar, |i?(0)|2 r\l-r)z'' 



dzd COS 6 \ yj >p g^ j^3 1^^ -^ (1 — rYz"^ 

Jo \ otsKIJ'h) J 

3 Results and Discussions 

It is understood that SCET is only valid at the large z region, while NRQCD should be fine 
in the small and medium z region. Therefore in order to obtain a formula which can describe 
the J /if) production in the whole kinematic region, one shall interpolate smoothly between 
the NRQCD and resummed SCET results. Here we propose an interpolating formula 

da\nt ,^ n'^O'NRQCD , C^O'resum /nri\ 



dzd COS Q dzd cos Q dzd cos Q 

Obviously the NRQCD contribution vanishes in the limit of z = 1, and only the resummed 
contribution survives, while at the small z^ the NRQCD contribution dominates. In addition, 
if one does not do any expansion and resummation in SCET, ciresum should be replaced by 
(Jnrqcd and hence Eq. (|^ will reproduce the NRQCD result. 
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Figure 2: The momentuin distribution (a) and angular coefficient a{z) (b) for e~^e~ -^ J/4'99 
process. The dashed, dotted and dot-dashed curves correspond to the NRQCD, tree-level 
SCET and resummed SCET calculations, respectively. The solid curve is for the interpolated 
resummed results. 

The differential cross section for the J /ipgg production are restricted by unitarity, parity, 
and angular momentum considerations. Its polar angle dependence can be parametrized into 
the form [3] 



da 



dzd cos 6 



S{z)[l + a{z)cosH], 



:2r 



where the angular coefficient a{z) is generally limited in the interval —l<a{z)<l. This 
general form has been conffimed directly by the calculations in the framework of NRQCD 
[inilSlIll- From Eqs. (fTH|l . (fTTj) and (fTT^ . it is easy to find that the tree-level and resummed 
SCET results also keep the form of Eq. (PT|) . and furthermore, have the same coefficient 
a{z). As a natural result, the interpolated resummed cross section in Eq. ^I^ follows the 
same behavior. 

In our numerical estimation we use y/s = 10.58 GeV and M = 2mc = 3.0 GeV. For 
simplicity, we also normalize the cross section to a dimensionless quantity by a factor R = 
il28/3)asM'cy'elM\RiO)\ys\ 

In Fig. 2(a), we show the momentum distribution of the process e"'"e~ — >■ J/ipgg. The 
dashed, dotted, dot-dashed and solid curves correspond to the NRQCD, tree-level SCET, 
resummed SCET calculations and the interpolated resummed result, respectively. The 
NRQCD result is taken from Ref. ^E], while the tree-level SCET, resummed SCET and 
interpolated resummed results are obtained by integrating over the polar angle cos 6 in Eqs. 
fITTj) . (fT^ and (PUJ) . As a cross-check of Eq. (fT7|) . one can see that the tree-level SCET result 
coincides with the NRQCD one in the limit z —>■ 1. Comparing to the NRQCD calculation, 
the interpolated resummed momentum distribution is suppressed significantly not only in 
the large z region but also in the medium z region. For example, at z = 0.9 the ratio of 
the interpolated resummed cross section and the NRQCD cross section is about 0.4, while 
at z = 0.5 the ratio is still 0.6 which is not quite close to unit. However, the large suppres- 
sion might be overestimated. This is because of the large discrepancy between the NRQCD 
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and tree- level SCET results. Although at z = 1, the NRQCD and tree- level SCET results 
are exactly the same, which is guaranteed by the matching procedure, the tree-level SCET 
spectrum deviates very quickly from the NRQCD one as z departs from one. For instance, 
at z = 0.9, the tree-level SCET cross section is only half of the NRQCD one. This indicates 
that the tree-level SCET calculation may not be a good expansion of the NRQCD calculation 
even in the large z region. The resummed SCET cross section is entirely based on the tree- 
level SCET calculation, as shown in Eqs. ()17p and (fT^ . and therefore the over-suppression 
occurs after interpolating between the NRQCD and the resummed SCET contributions. 

The discrepancy between the NRQCD and tree-level SCET results can be further in- 
vestigated by the J/?/^ angular distribution. In Fig. 2(b), we show the angular coefficient 
a defined in Eq. PT|) as the function of z. The dashed, dotted and solid curves are for 
the NRQCD, tree-level SCET and the interpolated resummed results, respectively. The re- 
summed SCET result has the same a(z) as that of the tree- level SCET. a in the NRQCD 
calculation is around zero in the region z < 0.85, while falls off rapidly as z > 0.85. At 
2; = 1, a is about —0.85. In contrast, a in the tree-level or resummed SCET result almost 
does not change with z. This behavior provides some hint about why the tree-level SCET 
result does not match with the NRQCD one very well at large z. It was known that at the 
end point z = 1, only the scalar component of the gluon-gluon system is allowed, which gives 
a = —0.85 [iZj. Apart from the end point, other spin components should be involved and 
might give dominant contributions which increase a fast to be around zero with decreasing z. 
The leading SCET expansion in the small parameter A ~ a/1 — z, which gives rise to a scalar 
operator for the gluon-gluon system (Eqs. (0) and (jHl)), cannot describe the contributions 
from other spin components. This implies that the power counting rules of SCET might 
break down due to some yet unknown reasons. One possibility is that part of the power 
suppressed contributions might be kinematically enhanced significantly. If this were true, 
one would have to match onto SCET to the next-to-leading order in A, obtain the power 
suppressed operators and their coefficients, and then perform the resummation procedure. 
However this complicated calculation goes far beyond the purpose of this paper. 

It is well known that the scale of J/ip is a little awkward for the application of NRQCD, 
thus one might worry whether the discrepancy between the NRQCD and tree-level SCET 
results is just an illusion. The key observation here is that, the only dimensionless parameter 
in this process is r = M^/s. Therefore the normalized cross section a/R, which is dimen- 
sionless, should only depend on r. That means even in a model world in which J/ip could 
be chosen to be very heavy (for example 30 GeV), the J/ip momentum spectrum would still 
be the same as that showed in Fig. 2 if r = M'^/s is taken to be fixed by increasing the cm. 
energy ^/s correspondingly. That is to say, even in a model world in which the application of 
NRQCD is guaranteed by very massive J/tp and the perturbative treatment of jet function 
is guaranteed by the larger cm. energy ^/s, the large discrepancy between the NRQCD and 
tree-level SCET calculations would still be there. 

As we have emphasized before the similarity between the radiative T decay T —>■ 'jgg and 
the inclusive J/ip production 7* -^ J/tpgg, one might naturally ask whether there is similar 
trouble for the former case. As shown in Ref. ^], for the photon momentum spectrum. 




Figure 3: The momentum distribution (a) and angular coefficient a{z) (b) for e"'"e~ — > 
^ ~^ 199 process. The dashed, dotted and dot-dashed curves correspond to the NRQCD, 
tree-level SCET and resummed SCET calculations, respectively. The solid curve is for the 
interpolated resummed results. 



there is no significant discrepancy between the NRQCD and tree-level SCET results for 
the radiative T decay. However Ref. [12] has not investigated the angular distribution of 
photons. Considering the process e+e" -^ T ^ 'ygg at CLEO, T is transversely polarized in 
the cm. frame. Accordingly the ultrasoft matrix element, which is proportional to v^v^ —9^^ 
in Ref. [13] {v is the T velocity), now should change to be proportional to —g^f^ (Eq. (JU) by 
using vacuum-saturation approximation. Performing an analogous calculation as what we 
have done in the last section, we obtain the differential T decay rate in the tree-level SCET 



dV 



trSCET 



T^dzd COS 9 8 



;z(l + cos^6'). 



(22) 



where z = 2E^/Mj; and 6 is the scattering angle between the momentum of the outgoing 
photon and the electron beamline in the cm. frame. Here Fq is a normalization constant. 
Eq. (^^ indicates that the angular coefficient a{z) defined in Eq. ()2H1 is always equal to 
unity at any z in the SCET calculation. 

The resummed SCET momentum distribution is the same as that in Ref. ^3]; while we 
choose the interpolation way as in Eq. (^01) in order to give an interpolated resummed result 
for both momentum and angular distribution. 

In Fig. 3, we show the momentum distributions and the coefficient a of the radiative 
T decay as the function of z in the NRQCD (the dashed line), the tree- level SCET (the 
dotted line), the resummed SCET (the dot-dashed line) and the interpolated resummed 
calculations. The NRQCD result is taken from Ref. ^Hl- It is clear that although the 
momentum distribution in the tree-level SCET calculation is close to that in NRQCD, the 
large difference oia{z) still exists in the end point region. Similar with the J/1P99 production, 
the gluon-gluon system here also leaves only scalar component at the end point z = 1 [19]. 
Therefore the NLO matching onto SCET might also play an important role in this process. 

In this paper, we studied the coUinear suppression effect in the process e+e~ — *> J /ip99 and 
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performed a leading power calculation in SCET. We obtained the decreasing J /%p spectrum in 
the endpoint region, which comes from the Sudakov logarithms suppression, and combine our 
SCET result with the NRQCD calculation. We then showed the momentum and angular 
distributions for J /if) in SCET and compared with the NRQCD results. Surprisingly, we 
found that, even before the resummation over large logarithms, there already exists a large 
discrepancy between the SCET and NRQCD results in the endpoint region of J /ip spectrum. 
A similar discrepancy is also found in the angular distribution of the radiative T decay. 
Therefore it should be highly interesting to have further investigations, for example, including 
the power suppressed contributions, on these processes. 

Acknowledgments 

We would like to thank Kaoru Hagiwara and Deshan Yang for helpful discussions. ZHL is 
supported by the Japan Society for the Promotion of Science (JSPS). 

References 

[1] G.T. Bodwin, E. Braaten, and G.R Lepage, Phys. Rev. D 51, 1125 (1995) [Erratum- 
ibid. D 55, 5853 (1997)]. 



[2] 
[3] 

[4] 
[5] 
[6] 
[7] 



[9 
[10 



CDF Collaboration, F. Abe et al, Phys. Rev. Lett. 79, 572 (1997); 79, 578 (1997). 

P. Cho and A.K. Leibovich, Phys. Rev. D 54, 6690 (1996); F. Yuan, C.F. Qiao, and K.T. 
Chao, Phys. Rev. D 56, 321 (1997); S. Baek, J. Lee, H.S. Song, and P. Ko, J. Korean 



Phys.i Soc. 33, 97 (1998); S. Baek, P. Ko, J. Lee, and H.S. Song, |hep-p h/9804455 
K. Hagiwara, E. Kou, Z.-H. Lin, C.-F. Qiao, and G.-H. Zhu, |hep-ph/0401246[ 

E. Braaten and Y.Q. Chen, Phys. Rev. Lett. 76, 730 (1993). 

BaBar Collaboration, B. Aubert et al, Phys. Rev. Lett. 87, 162002 (2001). 

Belle Collaboration, K. Abe et al, Phys. Rev. Lett. 88, 052001 (2002); ibid. 89, 142001 
(2002) 

P. Pakhlov, Measurement of double cc production, the talk given on 22 - 29 March 2003, 
Les Arcs, France. 

S. Fleming, A.K. Leibovich, and T. Mehen, Phys. Rev. D 68, 094011 (2003). 

F. Maltoni and A. PetreUi, Phys. Rev. D 59, 074006 (1999); I.Z. Rothstein and M.B. 
Wise, Phys. Lett. B 402, 346 (1997). 



10 



[11] C.W. Bauer et al, Phys. Rev. D 63, 014006 (2001); C.W. Bauer et al, Phys. Rev. D 
63, 114020 (2001); C.W. Bauer and I.W. Stewart, Phys. Lett. B 516, 134 (2001); C.W. 
Bauer et a/., Phys. Rev. D 65, 054022 (2002). 

[12] C.W. Bauer et al, Phys. Rev. D 64, 114014 (2001). 

[13] S. Fleming and A.K. Leibovich, Phys. Rev. Lett, 90, 032001 (2003); Phys. Rev. D67, 
074035 (2003). 

[14] X. Garcia i Tormo and J. Soto, Phys. Rev. D69, 114006 (2004). 

[15] CLEO Collaboration, B. Nemati et al, Phys. Rev. D 55, 5273 (1997). 

[16] W.Y. Keung, Phys. Rev. D 23, 2072 (1981); J.H. Kiihn and H. Schneider, Phys. Rev. D 
24, 2996 (1981); Z. Phys. C 11, 263 (1981); V.M. Driesen, J.H. Kiihn, and E. Mirkes, 
Phys. Rev. D 49, 3197 (1994). 

[17] S.J. Brodsky, A.S. Goldhaber, and J. Lee, Phys. Rev. Lett 91, 112001 (2003); S. Dulat, 
K. Hagiwara, and Z.-H. Lin, hep-ph/040223. 

[18] K. Koller and T. Walsh, Nucl. Phys. B 140, 449 (1978). 

[19] X.-G. He, H.-Y. Jin, and J.P. Ma, Phys. Rev. D 66, 074015 (2002). 



11 



